Possible indication to the QCD evolution of double parton distributions? 
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For the first time the process- independent parameter of double parton scattering, cr^ p , has been 
measured newly in the DO experiment at the three different resolution scales. If we interpret the 
measurement as a decrease of the effective cross section with a growth of the resolution scale it can 
indicate the QCD evolution of double parton distributions. 
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I. INTRODUCTION 

The presence of multiple parton interactions in high 
energy hadron collisions has been convincingly demon- 
strated by the AFS [J, UA2 and CDF [| Collab- 
orations using events with the four-jet final state and 
later by the CDF Collaboration [H using events with the 
7 + 3 jets final state. Recently the DO Collaboration has 
measured the process-independent parameter of double 
parton scattering, er°g P , using the 7 + 3 jets events at 
the three different resolution (energy) scales [f| provid- 
ing new and complementary information on the proton 
structure. The possibility of observing two separate hard 
collisions has been proposed since long Q,and from that 
has also developed in a number of works @, fl EMI ■ A 
brief review of the current situation and some progress 
in the modeling with account of correlated flavor, color, 
longitudinal and transverse momentum distributions can 
be found in Ref. [Toj . 

Multiple interactions require an ansatz for the struc- 
ture of the interacting hadrons, i.e. correlations between 
the constituent partons. As a simple ansatz, usually, the 
two-parton distributions (of the parton momentum frac- 
tion) are supposed to be the product of two single-parton 
distributions times a momentum conserving phase space 
factor. In recent papers [l4| it has been shown that this 
hypothesis is in some contradiction with the leading loga- 
rithm approximation of pcrturbative QCD (in the frame- 
work of which a parton model, as a matter of fact, was es- 
tablished in the quantum field theories [HI). Namely, the 
two-parton distribution functions being the product of 
two single distributions at some reference resolution scale 
become dynamically correlated at any different scale of a 
hard process. 

The main purpose of the present letter is to show that 
these QCD dynamical correlations result effectively in 
a dependence of the experimentally measured effective 
cross section cr^ p on the resolution scale unlike naive 
accepted expectations. 

The paper is organized as follows. Section [TT| is de- 
voted to what is known from perturbative QCD theory 
on the two-parton distribution functions. The possible 
manifestation of correlations induced by QCD evolution 
is discussed in Sec. IIIII We summarize and conclude in 
Sec. [TV] 



II. DOUBLE PARTON DISTRIBUTIONS IN 
THE LEADING LOGARITHM 
APPROXIMATION 

In order to introduce the denotations and to be clear 
let us recall that, for instance, the inclusive differential 
cross section for the four-jet production due to the si- 
multaneous interaction of two parton pairs within one pp 
collision, cr^p, can be given by Q 

(T dp = ^2 J a ^ 2 J^ i Dp(xi,x 3 )Dp(x2,x i )dxidx2dx 3 dx4 : , 
q/g 

.(1) 

where Uij stands for the two-jet production cross section. 
The dimensional phcnomcnological parameter a a s in the 
denominator is a factor characterizing a size of the ef- 
fective interaction region of the hadron (the factor 2 is 
introduced due to the identity of the two parton subpro- 
cesses). The relatively small value of (<T e fi)cDF measured 
by the CDF Collaboration Q with respect to the naive 
expectation was, in fact, considered 0, Q as evidence 
of nontrivial correlation effects between partons in the 
transverse space. But, apart from these correlations, the 
longitudinal momentum correlations can also exist and 
they were investigated in Ref. [h}. The factorization 
ansatz is just applied to the two-parton distributions en- 
tering Eq. ([1]): 

D p (xi,Xj,Q 2 ) = Dp(x l ,Q 2 )D p (x j ,Q 2 )(l-x l -x j ), (2) 

where D p (xi, Q 2 ) are the single quark/gluon momentum 
distributions at the scale Q 2 (determined by a hard pro- 
cess). 

However multi-parton distribution functions sat- 
isfy the generalized Gribov-Lipatov-Altarelli-Parisi- 
Dokshitzer (GLAPD) evolution equations derived for the 
first time in Refs [la, [l?} as well as single parton dis- 
tributions satisfy more known and cited GLAPD equa- 
tions [H, HU . Under certain initial conditions these gen- 
eralized equations lead to the solutions, which are iden- 
tical with the jet calculus rules proposed originally for 
multiparton fragmentation functions by Konishi-Ukawa- 
Veneziano [l9j and are in some contradiction with the 
factorization hypothesis ([2]). Here one should note that 
at the parton level this is a strict assertion within the 
leading logarithm approximation. 
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where b = (33 — 2ri/)/127r in QCD, g(fi 2 ) is the running 
coupling constant at the reference scale /i 2 , rt/ is the 
number of active flavors, Aqcrj is the dimensional QCD 
parameter, the GLAPD equations read [H, EH 



dDj(x,t) 
dt 



E 



(3) 



They describe the scaling violation of the parton distribu- 
tions Dj (x, t) inside a dressed quark or gluon (i, j = q/g). 

We will not write the kernels P explicitly and will not 
derive the generalized equations for two-parton distribu- 
tions D J i lj2 (xi,X2,t), representing the probability that 
in a dressed constituent i one finds two bare partons of 
types ji and j2 with the given l ong itudinal momentum 
fractions x\ and X2 (referring to [14T - [l8| for details). We 
note only that their solutions can be represented as the 
convolution of single distributions [1(1 EL This convolu- 
tion coincides with the jet calculus rules [19( as mentioned 
above and is the generalization of the well-known Gribov- 
Lipatov relation installed for single functions [Hj] (the 
distribution of bare partons inside a dressed constituent 
is identical to the distribution of dressed constituents in 
the fragmentation of a bare parton in the leading loga- 
rithm approximation) . The obtained solution shows also 
that the double distribution of partons is correlated in 
the leading logarithm approximation: 



Df^{x u x 2 ,t) ^Df {x u t)Df (x 2 ,t). 



(4) 



Of course, it is interesting to find out the phenomeno- 
logical issue of this parton level consideration. This can 
be done within the well-known factorization of soft and 
hard stages (physics of short and long distances). As a 
result, the equations ^ describe the evolution of parton 
distributions in a hadron (h) with t (Q 2 ), if one replaces 
the index i by index h only. However, the initial condi- 
tions for new equations at t = (Q 2 = /i 2 ) are unknown 
a priori and must be introduced phenomenologically or 
must be extracted from experiments or some models deal- 
ing with physics of long distances [at the parton level: 



D 3 Ax,t = Q) = SijSix-l); Dj lj2 ( Xl ,x2it = 0) 



0]. 



Nevertheless the solution of the generalized GLAPD evo- 
lution equations with a given initial condition may be 
written as before via the convolution of single distribu- 
tions [Til, [13 ■ This result shows that if the two-parton 
distributions are factorized at some scale /x 2 , then the 
evolution violates this factorization inevitably at any 
different scale (Q 2 ^ fi 2 ), apart from the violation due 
to the kinematic correlations induced by the momentum 
conservation, which is the analogue of the momentum 
conserving phase space factor in Eq. @. 



For a practical employment it is interesting to know 
the degree of this violation. Partialy this problem was 
investigated theoretically in Refs. Qj], H3] and for the 
two-particle correlations of fragmentation functions in 
Ref. [2l|. That technique is based on the Mellin trans- 
formation of distribution functions and the asymptotic 
behavior can be estimated. Namely, with the growth of 
t (Q 2 ) the correlation term becomes dominant for finite 
xi and X2 and thus the two-parton distribution func- 
tions "forget" the initial conditions unknown a priori and 
the perturbatively calculated correlations appear. 

The asymptotic prediction "teaches" us a tendency 
only and tells nothing about the values of xi,X2, t(Q 2 ) 
beginning from which the correlations are significant. 
Naturally numerical estimations can give an answer to 
this specific question using the CTEQ fit [22| for single 
distributions as an input. The nonperturbative initial 
conditions D ] h {x,Q) are specified in a parametrized form 
at a fixed low-energy scale Qo = fi = 1-3 GeV. The par- 
ticular function forms and the value of Qq are not cru- 
cial for the CTEQ global analysis at a flexible enough 
paramctrization. The results of numerical calculations 
were obtained in Ref. UM for the ratio: 



R(x,t) 



D" 

p(QCD,corr.) 



(X!,x 2 ,t) 



D 9 p (x 1 ,t)D 9 p (x2,t)(l-x 1 ~x 2 y 



(5) 



At the hard process scale of the CDF measurement [j] 
(Q ~ 6 GeV) the ratio ([5]) is nearly 10% and increases 
right up to 30% at much higher scale (Q ~ 100 GeV) 
for the longitudinal momentum fractions x < 0.1 acces- 
sible to these measurements. For the finite longitudinal 
momentum fractions x ~ 0.2 -r- 0.4 the correlations may 
increase right up to 90%. They become important for 
almost all x with growing t in accordance with the pre- 
dicted QCD asymptotic behavior [13, H^. The more com- 
prehensive analysis of double parton distributions based 
on the direct numerical integration of generalized equa- 
tions can be found in the quite recent paper [23j (see 
also Ref. [ll| for estimation therein of joint interaction 
probability due to evolution in comparison with the fac- 
torization contribution). 



III. POSSIBLE MANIFESTATION OF 
CORRELATIONS INDUCED BY QCD 
EVOLUTION 

As an effect of evolution, the double distribution func- 
tions become strongly correlated in longitudinal momen- 
tum fractions at large Q 2 and finite x as mentioned 
above. On the other hand, the indications from the ex- 
perimental observation of double scatterings at CDF [3| 
are not in favor of strong correlations effects in longi- 
tudinal momentum fractions. The most likely reason 
is that the studied kinematical domain, with relatively 
small x values and low resolution (energy) scale, is far 
from the asymptotic QCD prediction. The possibility of 
testing double collisions at much higher resolution scales 
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will open an opportunity of probing the correlations pre- 
dicted by the QCD evolution directly. For this purpose, 
the double parton cross sections of the equal sign W pair 
production (with a high resolution scale) in pp collisions 
at 1 TcV < y/s < 14 TeV are considered in Ref. [H[ . As 
a main result, the contribution of the term with corre- 
lations in the equal sign W pairs production might con- 
tribute almost 40% to the cross section at yfs = 1 TeV 
and about 20% at 14 TcV. 

Here we would like to bring attention to the DO mea- 
surement of the process- independent parameter of double 
parton scattering, a^ p , done as a function of the second 

(ordered in the transverse momentum px) jet pr, P^ t2 , 
that can serve as a resolution scale. It is shown in Fig. Q] 
(see also Fig. 11 from Ref. @) and can be considered as 
a first inderect manifestation of the evolution effect since 
<7^ p shows a tendency to be dependent on the resolution 
scale. At first glance the dimensional parameter a e g en- 
tering into Eq. (fTJ) contains the information related with 
the non-perturbative structure of the proton and corre- 
sponds to the overlap of the matter distributions in the 
colliding hadrons and therefore it is independent of the 
hard process scale a priori. However the experimental 
effective cross section tr°g P is not measured directly but 
is calculated (extracted) using the normalization to the 
product of two single cross sections: 



a 



7+3.7 
DPS 



exp-i-1 



(0) 



in both the CDF and DO experiments. Here a 10 and a 33 
are the inclusive 7+ jet and dijets cross sections, crj^p^ 
is the inclusive cross section of the 7 + 3 jets events pro- 
duced in the double parton process. The factor 2 before 
cr^ p is not needed in this case of distinguishable scat- 
terings unlike the four-jet process. It is worth noticing 
that the CDF and DO Collaborations extract a e ^ with- 
out any theoretical predictions on the 7+ jet and dijets 
cross sections, by comparing the number of observed dou- 
ble parton 7 + 3 jets events in one hard pp collision to 
the number of 7 + 3 jets events with hard interactions 
occurring in two separate pp collisions. 

At such a normalization ((6]) and the presence of corre- 
lations in the two-parton distributions the experimentally 
extracted dimensional factor cr°^ p will be different from 
one cr e ff incoming in Eq. (TTJ). Indeed, in accordance with 
QCD evolution instead of the factorization ansatz © we 
can write [3 Hi]: 

D t p 3 {x l ,x 3l t) = D p { Xi ,t)D 3 p {x h t){l -an- Xj)[l + 




FIG. 1: Effective cross section a^ p measured in the three p^ t2 
bins at the DO experiment [||. The solid (k = 0.5) and dashed 
(k = 0.1) lines are the results from Eq. (|ll|l at p^ = 22.5 
GeV and = 16.3 mb. 



Dp(xi, 0)D p (xj, 0)(1 — Xi — Xj) which is supposed to be 



factorized at the reference scale t = (Q 



and 



Dp(xi,t)D 3 p (x 3 ,t)(l - Xi ~ Xj)Rz (xi,Xj,t) (9) 



is a particular solution of the complete equation with 
zero initial condition R% (xj,Xj,0) = 0. At t > this 
solution is always positive as it is the integral convolu- 
tion of positive single distributions with the kernels of 
nonhomogeneous part of the evolution equation. These 
kernels (probabilities) arc defined without negative 5- 
function rcgularization and therefore are always positive 
unlike the kernels of homogeneous part which have neg- 
ative contributions. In the kinematical range of interest 
for the actual case (we never exceed x = 0.1 practically) 
the contribution of the term (xi , Xj , t) to Eq. ([7]) is 
negligible [HI]. Substituting Eq. in Eq. (UJ, and ne- 
glecting the factor (1 — Xj — Xj) as it is usually done at 
relatively small Xi, Xj, and using Eq. (J6j) we derive 



[aeffJ-^l + A^)], 



(10) 



where A(i) is a some positive contribution induced by 
the correlation term R£ (xj, x 3 ■, t). This contribution in- 
creases OH] with a growth of the resolution scale t(Q 2 ) 
and can be right up 40% at Q = M w = 80.4 GcV fol- 



R\ J (xi,Xj,t) + R%(xi,Xj,t)], (7) lowing Ref. [24|. Taking into account that namely cr off 



where 



D p {xi,t)D 3 p {x J ,t)(l -x t - Xj )[l + R^ix^x^t)} (8) 

is the solution of the homogeneous generalized GLAPD 
evolution equation with the given initial condition 



is the true process-independent dimensional parameter 
we conclude: the experimentally extracted effective cross 
section must decrease with the growth of the reso- 
lution scale. The experimental data in Fig. [1] show just 
such a tendency in spite of large experimental uncertain- 
ties which do not allow the DO Collaboration to make 
this conclusion. 
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Unfortunately, the calculation of correlation contri- 
bution A(t) to the effective cross section tr°g P is not 
yet possible in the framework of some kind of existing 
event generators. For instance, the double parton scat- 
terings are implemented in the Monte-Carlo generator 
PYTHIA 0, [27J taking into account some correlations 
which, however, are not quite adequate for the case of 
evolution effects under consideration (unlike the theoret- 
ical investigations in Ref. [ll|). The implementation of 
the QCD evolution of two-parton distribution functions 
in some Monte-Carlo generator, as this was done for sin- 
gle distributions, is not a trivial task. However we can 
"predict" a functional form of p^ t2 -dependence of cr^ p : 

a^^^ + kH^/rto}- 1 (11) 

inspired by the explicit expression for the correlation 
term [l4[ and the evolution variable t. The results from 
Eq. (fTTj) are also shown in Fig. 1 at k = 0.1 (dashed 
line) and k = 0.5 (solid line) to illustrate the two possi- 
ble slopes of cr°g P in the observation region: visually im- 
perceptible (practically constant) and noticeable falling 
slopes. The normalization point p^q 2 = 22.5 GeV with 
<7° ff = 16.3 mb is fixed to reproduce the experimental 
value in the central p!^ t2 bin. 

IV. CONCLUSIONS 

We argue that the QCD dynamical correlations result 
effectively in the dependence of the experimentally ex- 
tracted cr°^ p on the resolution scale unlike naive accepted 
expectations. The measurements covering a larger range 
of the resolution scale variation with a smaller uncer- 



tainty are needed to see the evolution effect more dis- 
tinctly. 

In order to investigate the more delicate characteris- 
tics of double parton scatterings (distributions over var- 
ious kinematic variables with various kinematic cuts) it 
is also desirable to implement the QCD evolution of the 
two-parton distribution functions in some Monte Carlo 
event generator as this was done for the single distri- 
butions, for instance, within PYTHIA [26|. It is worth 
noticing once more that the evolution of the two-parton 
distribution functions has the same confidence status as 
the well-established evolution of the single distribution 
functions in the framework of the leading logarithm ap- 
proximation of perturbativc QCD. Therefore the exper- 
imental direct or indirect observation of this evolution 
effect is significant to answer to many challenging ques- 
tions of yet poorly-understood aspects of QCD. 
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